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Abstract 

Recent experimental data for the complete wetting behavior of pure ^He and of ^He-^He mixtures 
exposed to solid substrates show that there is a change of the corresponding film thicknesses L 
upon approaching thermodynamically the A-transition and the tricritical end point, respectively, 
which can be attributed to critical Casimir forces fc- We calculate the scaling functions of fc 
within models representing the corresponding universality classes. For the mixtures our analysis 
provides an understanding of the rich behavior of "d deduced from the experimental data and 
predicts the crossover behavior between the tricritical point and the A-transition of pure ^He 
which are connected by a line of critical points. The formation of a 'soft-mode' phase within the 
wetting films gives rise to a pronounced maximum of fc below the tricritical point as observed 
experimentally. Near the tricritical point we find logarithmic corrections ~ L~'^(lnL)^/^ for the 
leading behavior of dominating the contributions from the background dispersion forces. 

PACS numbers: 05.50.+q, 64.60.Cn, 64.60.Kw, 67.40.Kh 
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I. INTRODUCTION 



There is growing experimental evidence for the analogue of the electromagnetic Casimir 
effect [l| in various critical condensed matter systems {3, 3, 4, 5, 6, 7|. In wetting experiments 
the confinement of critical fluctuations within an adsorbed liquid film gives rise to an effective 
Casimir force fc between the substrate-liquid and the liquid-vapor interfaces of the liquid 
film j^, y, [l^. Near the critical end point of the liquid the emerging Casimir force adds to 
the omnipresent dispersion forces and thus leads to a change of the thickness of the complete 
wetting film. From this response one can infer the Casimir force by subtracting the effect of 

;he critical end point with temperature 



111 ] this force fc per unit area and in 



the background forces which varies smoothly near 
Tc. In accordance with finite-size scaling theory 
units of ksTc can be expressed in terms of a universal scaling function i); its shape depends 
sensitively on the type of boundary conditions (BC) [9] and thus on the surface universality 



classes the confining surfaces belong to 



121. 



Capacitance measurements of the equilibrium thickness of ^He wetting films near the 
superfiuid temperature Tx of the critical end point of the A-line {2, 7| quantitatively support 
the theoretical predictions of fc for the bulk universality class of the XY model with sym- 
metric Dirichlet-Dirichlet BC (O, O) forming the so-called ordinary (O) surface universality 
class Q]. Such BC correspond to the case that the quantum-mechanical wave function of 
the superfiuid state vanishes at both interfaces, giving rise to an attractive Casimir force 
{fc < 0) 0, Q. However, the available theoretical results have a limited range of apphcabil- 
ity, i.e., T > Tx and T <^Tx- Above and at Tx explicit field-theoretical calculations within 



the e-expansion scheme are available 



13 



14| . For temperatures well below Tx there are cal- 



culations which take into account capillary-wavelike surface fluctuations in the asymptotic 
limit of thick films, predicting a levelling off of the scaling function for large negative scaling 
variables [l^, i.e., T -C Tx, in qualitative agreement with the experimental observations. So 
far there are no theoretical results available for the critical region below Tx which provide 
an understanding of the deep minimum of the experimental scaling function (ca. 20 times 
deeper than its value at Tx). 



^He-'^He mixtures near their tricritical end point (see Fig. 12 in Ref. 14 [ 



are another 
, l5| . The tri- 



critical system for which wetting experiments have been performed recently 
critical end point with temperature Tt is the point in the ^He-'^He phase diagram where the 
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line signalling the onset of superfluidity joins the top of the two-phase coexistence region 
for phase separation into a ^He-rich superfluid phase and a '^He-rich normal phase. The 
mixture belongs to a bulk universality class different from that one of pure ^He and, because 
its upper critical spatial dimension d* equals 3, the actual physical system is character- 
ized by rational mean-field critical exponents (up to logarithmic corrections) ly, ll7|. The 



capacitance measurements of the wetting film thickness of the mixture reveal a repulsive 
Casimir force fc around the tricritical end point which suggests non-symmetric BC for the 
superfluid order parameter (OP). The probable physical mechanism behind such a BC is 
that within ^He-^He wetting films a ^He-rich layer forms near the substrate-liquid interface, 
which may become superfluid already above the line of onset of superfluidity in the bulk 
[isl whereas the lighter ^He has a preference for the liquid-vapor interface. Thus the two 
interfaces impose a nontrivial concentration profile which in turn couples to the superfluid 
OP. 

For this system, recently \^\ we briefly reported explicit model calculations which demon- 
strate that the concentration profile indeed induces indirectly non- symmetric BC for the su- 
perfluid OP. For symmetry-breaking (+) BC at the substrate-liquid interface and Dirichlet 
(O) BC at the liquid- vapor interface we calculated the Casimir force and found a semiquan- 
titative agreement with the experimental data given in Ref. [J]. Moreover, we formulated 
theoretical predictions for the behavior of fc in the crossover regime between the tricritical 
point and the A-transition of pure ^He which are connected by a line of critical points and 
provided the universal leading behavior of the Casimir force at the tricritical point. 

The purpose of the present study is to elucidate the details of the two complementary 
approaches used in Ref. 19| and to extend them in order to obtain new results both for 
the tricritical ^He-^He mixture and the critical pure ^He. The presentation is organized 
as follows: In Sec. [Til we discuss the universal properties of the Casimir force. As already 
mentioned above, for the present tricritical behavior the upper critical dimension d* equals 
3 and therefore the thermodynamic functions of three-dimensional systems exhibit power- 
law behaviors with critical exponents taking their classical values. However, logarithmic 
corrections to the mean-field (MF) behavior are expected under experimental conditions 



171]. Using field-theoretical methods and renormalization-group (RG) analyses we obtain 



the leading asymptotic behavior of the Casimir force at the tricritical point. As a function 
of the film thickness L it has the form of a power law multiplied by a fractional power 
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of a logarithm and by the universal Casimir amphtude. In addition, we also derive the 
form of the finite-size scaling for the Casimir force in the vicinity of the tricritical point. 
As expected also the arguments of the associate scaling function acquire logarithmic 
corrections. These scaling functions are compared with the ones deduced from the experi- 
mental data in Ref. [^]. In Sec. IIIII we study within mean- field theory (MFT) films of the 
lattice vectoralized Blume-Emery-Griffiths (VBEG) model which belongs to the same 
universality class as the ^He-'^He mixture but is simple enough to allow for systematic studies 
of fc along all thermodynamic paths followed in the wetting experiments of Ref. ^. This 
facilitates the exploration of the crossover between the tricritical point Tt and the line of 
critical points and the coexistence region below Tf. This enables us to follow the Casimir 
force upon continuously switching the bulk universality class (from tricritical to critical) by 
changing the concentration of the ^He-^He mixture. The scaling functions corresponding 
to thermodynamic paths of constant concentration of the two components of the '^He-^He 
mixtures are calculated and compared with the corresponding experimental data in Ref. {4 1 . 
As a limiting case the VBEG model can describe also a film of pure ^He which is studied in 
Sec. II VI within MFT. The scaling function of the corresponding Casimir force is obtained in 
the critical region below Tx and compared with that one extracted from the experimental 
data in Ref. |2|. We also compare these results with the mean-field predictions which fol- 
low from the Landau- Ginzburg theory in the film geometry with suitable BC. In Sec. iVj we 
discuss the theoretical results obtained within the VBEG model and assess their relevance 
for interpreting the experimental data. We conclude with a summary and an outlook in 
Sec.lVD 



II. UNIVERSAL PROPERTIES 



For film geometries, in this section we investigate the universal properties of the Casimir 
force near tricriticality. In general two-component systems are characterized by the ordering 
density $ and its conjugate field h, and by a non-ordering density x and its conjugate 
field A. For liquid ^He-^He mixtures, x, and A correspond to the superfluid OP, to 
the ^He concentration and to the difference between the chemical potentials of the ^He 
and ^He components, respectively, whereas the field h conjugate to the superfluid OP is 
experimentally not accessible. 
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A. Scaling function from Landau- Ginzburg theory 



In order to capture universal properties we consider the standard dimensionless 0{n)- 
symmetric Landau- Ginzburg (LG) Hamiltonian for a tricritical system in the film geometry: 



where L is the film thickness, $ is the n-component order parameter OP (n = 2 corresponds 
to the XY universality class), and z is the coordinate normal to the confining surfaces; tq, 
Uq, and Vq are bare coupling constants depending, inter alia, on the temperature T and the 
non-ordering field A. ro(uo) = and mq > define the critical line, whereas at the tricritical 
point one has ro = uq = 0,vo > 0. The semi-infinite version of Eq. ([T]) has been studied in 
the context of surface critical behavior [21]. In the film geometry the Casimir force per area 
A of the cross section of the film and in units of /cbTj, 

fc^-{dr/dL) = (2) 



is given by the thermal average of the stress tensor component Tzz jo]: 

r{L) = (/ - m/iksz) (3) 

where / is the total free energy of the film per volume V = LA and fb is the bulk free energy 
density. For large L the excess free energy can be decomposed into surface and finite-size 
contributions: f^^iL) = fg^i + fs,2 + ^f{L). The stress tensor is given by |9| 

%^ = d,^ ■ d,^ - 6,jC -{d- 2)/{A{d - l)Mdj - 5.,V2)$^ (4) 

where C is the integrand in Eq. ([1]). In what follows we assume $ = 0, . . . , 0), i.e., we 

neglect helicity. For non-symmetric BC its relevance for the behavior of the Casimir force is 
not clear because the OP has the additional freedom to rotate across the film by a position 
dependent angle (f){z); the analyses of the role of helicity is left for future research. Within 
MET for the LG Hamiltonian, the determination of the tricritical Casimir force in the film 
geometry starts from the Euler-Lagrange equation 

As discussed in Sec. HI (+, O) boundary conditions, with the substrate at 2; < and vapor 

at z > L, 

m(0) = +00 and m(L) = (6) 



are supposed to mimic the experimental system of ^He-^He wetting films as studied in 
Ref. 4|. According to Eq. (jlj) the stress tensor component 7^^ evaluated within MFT and 
with $ = {m{z), 0, . . . , 0) for the OP (in the present MF approach we omit the brackets (•) 
indicating the thermal average) yields 

r.. = l{m'{L)r. (7) 

In deriving this expression we have used the property that %z = const throughout the film 
including the surfaces and we have chosen zq = L a.s the point of reference at which 7^^ is 
evaluated. Accordingly, the first integral of Eq. is given by 

{m\^oi^)r = 2%, + romloi^) + ^mt,o(^) + ^<oi^)- (8) 

Dimensional analysis yields that, at the upper critical dimension d = d* = 3, 
■m{z, L, To, Uq, Vq) can be expressed in terms of a dimensionless scaling function ^p+^o'- 

m+^oiz,L,ro,Uo,VQ) = (^-^^ L'^^^^ip+^oiz/ L,roL^ ,UoL;vo), (9) 

where Vq is dimensionless. Similarly, within this approach the normalized Casimir force can 
be expressed in terms of a dimensionless scaling function 'i9+,o: 

— 1 /2 

%z = fc{L,ro,Uo,Vo) = L-'^l'^^{roL\uoL,vo). (10) 

Equation ([8]) can be written in terms of these scaling functions v^+^o and ^+,o- 

/ 5 \ l/2 

i'P'-,,oi^)r = <o + roL' ^loi^) + ^—j uoL y,t,o(^) + vl,oi^), (H) 
where x = z/L. In turn, Eq. ( |TT|) can be integrated directly yielding the implicit equation 

(12) 

/ \ 1/2 







for the scaling function i)^Q{rQL'^,uoL,VQ). Note that the coupling constant vq > re- 
mains undetermined within mean-field theory and enters into "t^^o only in the combination 
Vq uqL. Under the renormalization group flow, at the upper critical dimension {d* = 3) 
the renormalized coupling constant v associated with vq tends to its fixed point value v* = 0. 
This RG flow generates logarithmic corrections to scaling due to the singular dependence of 
the scaled quantities on v (see, e.g., Eqs. Qj and ffTOl) ). With the transformation 

^=ra''% (13) 
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for the integration variable one can rewrite Eq. (fT2l) in the more convenient but still implicit 
form 

^ ^'""^ Jo y/l + ap^ + bp^+p<^' 
where the dimensionless parameters a and b are given by 

a = roL^^Z)-'^' and ^ = (^) «o)"'^' • (15) 

The numerical evaluation of the scaling function amounts to the following steps: (1) spec- 
ifying values for a and b, (2) evaluating "^^q fro^i Eq. (IT^ . (3) determining the values of 
the two scaling variables roL"^ and ^^"^uqL from Eq. ( JTSl) . 

From the symmetry properties of the order-parameter profile for the symmetry breaking 
opposing boundary conditions (+, — ) it is obvious that within MET the force for a film of 
thickness L in this case can be obtained from Eqs. f|T^ . holding for (+, O) BC, and ffTOl) by 
replacing L L/2 therein. This implies 'd^^{x,y) = 8^9^^q(x/4, ?//2). In the following we 
shall refer only to the (+, O) BC and drop the corresponding index. 

The precise dependence of tq and uq on the thermodynamic fields T and A is not known. 
Therefore it is not obvious how to follow in terms of these variables a specified path in the 
phase diagram such as the experimental path of fixed '^He concentration. However, assuming 
that tq and Uq are analytic functions of T and A in the neighborhood of the phase transition 
one can use the expansion [17j: 



= A(A)(T-T,(A)) + 0((T-T;,(A))2) and = i?(A) + 0((T - r,(A))), (16) 



where Tx(A) denotes the critical temperatures of the line of continuous phase transitions as 
a function of A, and B{A) and ^(A) are positive and non-zero on this line; -B(A) = at 
the tricritical point. 

In view of comparisons with experimental data, which we shall discuss later, it is useful to 
mention the relation between the parameters ro and uq and the experimentally controllable 
thermodynamic fields T — Tt and A — A^ where A = A^ at the tricritical point and Tx{At) = 
Tt. These "deviating fields" are not the proper scaling fields and it was shown (22! that a 
suitable (dimensionless) choice is provided by 

t = {T-T,)/T, and ^ = (A - AO/(M) + a't, (17) 
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where a' is the slope of the hne tangential to the phase boundary at the tricritical point. 
Thus for t — > with (7 = the tricritical point is approached tangentially to the phase 
boundary. Instead of t one could also use a scaling variable which is orthogonal to the 
loci (7 = 0; this would not affect the leading singular behavior for t,g ^ [17]. Near the 
tricritical point B{A), A{A), and Tx{A) can be expanded in terms of g and t. Using Eq. fll7p 
one has T -Tx{A) = T -Tt + {a'kB)-\A - At) + 0((A - Atf) = {a'Y^Ttg + 0((A - Atf) . 
Expressing A and T as a function of t and g one finds: 

To = A^g + A^t' + Oig^gt) and = B^t + B^g + 0{gt, g^^) (18) 

where Ai > 0,Bi > 0, A2, and B2 are constants. Due to the analytic structure of Eq. (fT6|) 
and because (A — A^) = kBTt{g — a't) the coefficient tq does not contain a term linear in t 
so that Mo ~ t + O(t^) if tq = 0. On the other hand ~ (? + 0{g^) if uq = 0. 



B. Logarithmic corrections at T = Tt 



At the tricritical point a = 6 = Eq. reduces to 

/ dp/Vl +p6 ~ 1.40218. 
Jo 



(19) 



Accordingly, in units of AksTt, the MET result for the tricritical Casimir force in the 
case of (+, O) BC is (see Eq. 1^) 



/*7~ 2.75684 



1/2 



(20) 



In c? = 3 — e the MET result at tricriticality (Eq. (!20l) ) yields the leading contribution in a 
perturbation series, i.e.. 



90 



(21) 



where both tzz = 2.75684L~^ and {%z)i do not depend on vq. After removing ultraviolet 
singularities via renormalization (R) the asymptotic scaling behavior of fc follows from 
substituting the renormalized v by the appropriate fixed-point value f * oc e. At d = d*, and 
under spatial rescaling by a dimensionless factor i, v fiows to its RG fixed point value v* = 



according to 



2l| 



2407r^ 
3n + 22 



1 

In^ 



+ c- 



In I In J 

h?i 



+ 



(22) 
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where v{i) is the running couphng constant with the initial condition vii{i = 1) = vr. With 
the rescahng factor i = Iq/L, where lo is a microscopic length scale of the order of a few 
A, this yields a logarithmic correction to the power-law dependence on L of the tricritical 
Casimir force: 



fc ^ 0.54(3n + 22)i/2^1n(L//o))'/'l>- 



c ln|ln(L//o)| 
2 |ln(L//o)| 



(23) 



Determining the constant c requires to extend the analysis in Ref. [21[ which is left for future 
research. Gaussian fluctuations give contributions of at least 0{v^) which are therefore of 
order and thus subdominant (see Eq. ( l23i) ). We compare Eq. ( 123|) for n = 2 with the 
data obtained by Garcia and Chan 4] for their experimental value of L 520A and for 
Iq ^ 1.3 A, the experimental value of the correlation length amplitude Cn = C(^)/N with 



23|. For these values 



z/f = 1 for concentration fluctuations below Tf in the superfluid phase 
Eq. ([23]) predicts 

i9t = fcL^ ^ 6.96 (24) 

whereas = 8.4 ± 1.7. The value of the theoretical function i)t at Tj, with Zq between 1 
and 2 A, is in reasonable agreement with the measured ^9^^^. In order to extract the actual 
value of the universal Casimir amplitude (i.e., the numerical prefactor 0.54v^28 = 2.86 in 
Eq. ( l23l) ) the experimental data call for a re- analysis based on the functional form given 
by Eq. ( l23i) . which renders the comparison independent of the choice for Iq, and requires to 
take into account the correction terms given in Eq. fl23|) . We want to emphasize that the 
tricritical Casimir force offers the opportunity to observe the so far experimentally elusive 
logarithmic corrections associated with tricritical phenomena. We note, that at tricriticality 
the Casimir force fc{L oo) dominates over the background dispersion forces. This differs 
from the case of critical Casimir forces for which both contributions decay with the same 
power law. It is interesting that the Casimir amplitude for the present (+, O) BC is the 
same as for (+, +) BC considered in Ref. 24l |. 



C. Logarithmic corrections to the scahng function 

The scaling properties of the Casimir force follow from the renormalized finite-size con- 
tribution to the excess free energy (Eq. ([3])). For carrying out the renormalization procedure 
of this quantity two aspects are relevant. First, for the film geometry, the width L of the 
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system is not renormalized ll|. Second, in the renormalized (R) finite-size contribution 
to the free energy Sf{L) (see the text before Eq. (jl])) the contributions from the additive 



counter terms cancel and one has 



5f^{r, u, v; fx, L) = 5f{ro, Uq, Vq; L) 



12 



(25) 



where the bare quantities Mq, tq, and Vq are expressed in terms of renormahzed ones r, m, and 
f ; /i is an arbitrary momentum scale. Since we are not considering correlation functions at 



the surface, all renormalization factors Z are the same as those in the bulk 
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2l|: 



ro = ZrV + U^/i ^"P, Uq = ZyU, = ^TX^Z^V^ 



(26) 



where the dimensions of the coupling constant are [ro] = yU^, [uo] = 1^!^^'' and [t>o] = 
Explicit perturbative results for the tricritical bulk renormalization functions Z^., P, and 



Z^ are known (see, e.g., Refs. 17|, l2l|). From Eq. fl25|) the RG equation can be derived in 
a standard fashion by exploiting the fact that 5/(ro, uq,vq\ L) is independent of fi. Because 
in Eq. ( l25i) there are no additive renormalization terms it follows that 6f^{L) satisfies the 
following homogeneous RG equation [l2|: 



K=r,u,v 



sni) = 



(27) 



where I3^{r,u^v\t) = fid^lon and 9^|o denotes derivatives with respect to /i at fixed bare 
interaction constants for k = r,u,v. The RG equation is solved by using the method of 
characteristics (see, e.g., Ref. |26|]): 



6f^{r\ u, V, /x; L) = 5/" (f'(£), u(£), v{i); fii; L) 



(28) 



where ^ is again a dimensionless spatial rescaling factor, k(£) are the running coupling 
constants with the initial condition R{1) = k, and due to the form of the renormalization of 
To (see Eq. fl26|) ) the new variable r' is given by 



17 



2l| 



r' = r + w{v, {j,)u^ 



(29) 



0^ 



For an explicit expression of w{v,fi) see Refs. [l7|, l2l|]. Equation fl28l) summarizes the RG 
transformation and the non-renormalization of L. Using dimensional analysis one obtains 

f'{i) u{i) v{i) 



6f''{r',u,v,fi-L) = {f,iY'-''>6f^ 



l,Lfii 



(30) 
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The desired asymptotic scaling behavior of Sf^ follows by substituting on the rhs of 
Eq. fl30|) the appropriate fixed-point values for the running couplin g co nstants f', u, and v. 



The infrared stable fixed point lies at v* = (240/ (Sri + 22) )e + O(e^) 2l|]. Upon approaching 
the upper critical dimension f * — ^> and for e ^ the relevant logarithmic corrections to 
the classical exponents are generated by the flow of the coupling constants under the RG 
transformation i ^ 0. In the limit £ — 0, v{i) is given by Eq. (122|) . The running variables 
f'{l) a.nd u(l ) ca n be written as f'{i) = Er{i]v)r' and u{i) = Eu{i;v)u. A straightforward 



analysis 
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2l| shows that Eri£;v) const and ^„(£;f) ~ | ln£|-2("+4)/(3n+22) £ ^ q. 
Choosing /i = I/Iq, fiiL = i{L/lo) = 1, and omitting the constant factor Er we obtain the 
following scaling form for 6f: 

6f^{r',u,v,fi- L) = L-Hf\r'L\uL\ HL/h)\-^^'^^'^'^^^+'^\ \ ln(L//o)|-^ 1, 1). (31) 

Due to Eq. ([2]) the scaling form for the Casimir force follows from Eq. fl3T|) as: 

fc{r\u,v- L) ^ L-''e{r'L\uL\ HL/l,)\-^^-+'^/^'^'^+^^\ \ HL/l,)\~'). (32) 

The scaling function 9 is given in terms of 5/^(2:1,^2,-23,1) as 9 = 25 + 
2zi{d6f^/dzi) — Z2{d5f^ / 8x2). The higher-order terms neglected in Eq. (l32l) are of the 
form L-3(ln(L//o))-i(2(n + 4)/(3n + 22))z2{d5 f"^ / 8x2) + L'^hi^L/Qy^z^^ddf'^/dz^) + 
L-3(-l + 2c(ln I ln(L//o)|)/(ln^ \{L/lQ)\){d5f^ /dz^). The third term in the latter expression 
stems from the correction to 2:3 (see Eqs. ( I30l) and ( l22l) ). 

At the upper critical dimension the asymptotic critical behavior obtained from the pertur- 
bative RG calculations within the Gaussian approximation is expected to be exact. However, 
at the lowest order, often referred to as renormalized mean-field theory (RMF) - which yields 
the free energy correctly with the leading logarithms - one neglects the contributions stem- 
ming from the Gaussian fluctuations and replaces the scaling function by its mean-field-like 
form but with the rescaled arguments. 

Applying this reasoning to the free energy we use the mean-field result given by Eqs. (fTOj) 
and f|T2l) with tq replaced in favor of r' according to Eq. fl29|) with w{v{i), fi{i)) const 
as £ ^ 0, Mo replaced by m| ln(L//o)r^("+^^/(^"+^^\ and vq replaced by ((2407r2)/(3n + 
22)) I ln(L//o)|^"'^ to obtain at lowest order: 

1 /2 

(33) 
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In the following we want to compare the behavior of the MF and RMF expression for the 
Casimir force. As we have already stressed before, fc calculated within the MF approach 
depends on the non- universal and dimensionless parameter vq (see Eq. fllOl) ). Upon compar- 
ing with experimental data this parameter can be used to fit the amplitude of the Casimir 
force, because Vq appears (albeit not exclusively) as a prefactor of the scaling function. 
The factor Vq , which multiplies the coupling constant uq (see the text after Eq. (fT2|) ). is 
absorbed in the definition of the scaling variable. 

In Fig. [U we have plotted two curves: (1) d^^iroL"^ = 0,y^^^) = fcL^ivo/90Y/^ as a 
function of y^^^ = {5/{2vo)Y^'^uoL (see Eq. ( |T2|) ). Here, the non-universal factor f o is 
absorbed in the definitions of the scaling function and of the scaling variable. As already 
mentioned before ~ t if tq = 0, so that UqL ~ tL. (2) fcL^ = {}Rmf (^^^yRUF^ ^ 
(28/(87rV3))i/2(ln(L//o))'/'t9^^^(0,?/^*^^) (for n = 2), where y^^^^ = uLihi{L/k)f'^\ 
Here, renormalization fixes the amplitude of the Casimir force replacing the non-universal 
prefactor (|^) (see Eq. ffTOj) ) of the scaling function by the amplitude and the logarithmic 
correction to the L dependence. The scaling variable y^^^^ includes the logarithmic correc- 
tion I ln(L//o)|~^^"'*'^''''^^"''^^^'' to u and the additional logarithmic term | ln(L//o)|^^^ stemming 
from the factor (5 / {2vo)y^'^ (see y^^^ and Eq. f|T2|) ). The numerical factor 7/(247r^) has been 
included into the definition of u. For comparison with experimental data this factor can be 
combined with the non-universal constant of proportionality between u and t. For the plot 
we have chosen the experimental value for L/Iq, i.e., 520 A/1.3 A. The shapes of both scaling 
functions are similar but the RMF result gives the correct value for the Casimir amplitude 
and the correct L-dependence of the scaling function. This should be helpful for interpreting 
experimental data obtained for different film thicknesses. 

In Fig. [2] we show the corresponding results for u = so that r' = r (see Eq. fl29l) ). r t, 
and rL^ ~ gL"^. We find that for u = both scaling functions decay much faster to zero 
than for r = 0. 

For To = one has m ~ t so that, up to the logarithmic corrections, the scaling function 
^MFj^g^ ^A/F-j gj^ould Correspond to the experimental curve d{tL) in Fig. [3] for the tricritical 
concentration ^. (We note that the argument of the experimental curve is given in units 
of A.) The solid line in this figure represents ^'^^{0,y^^^) suitably adjusted with respect to 



the parameter t>o such 
the experimental ones 



;hat the Casimir amplitude and the position of the maximum equal 

3- 
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III. VECTORALIZED BLUME-EMERY-GRIFFITHS MODEL. 



Based on the motivation provided in the Introduction, in this section we extend the 
VBEG model to the film geometry and study ^He-^He mixtures. 

A. The model 

We consider a three-dimensional slab of a simple cubic lattice consisting of L parallel 
(100) lattice layers with lattice spacing a so that L = La. Each layer has A = A/ a? sites, 
labeled z, j, . . ., which are associated with an occupation variable = 0, 1 and a phase Oi 

< 9i < 27r) which mimics the phase of the ^He wave function and thus renders the XY 
bulk universality class (ra = 2). A ^He (^He) atom at site i corresponds to = 0(1) so that 
in the bulk X = 1 — (ti) is the ^He concentration. Unoccupied sites are not allowed so that 
the model does not exhibit a vapor phase. Accordingly this model does not allow for the 
occurrence of a tricritical end point. However, we expect that the universal properties we are 
interested in are the same for tricritical points and tricritical end points. The Hamiltonian 
consists of bulk and surface contributions H = Hb + Hs with 

nb = -JY^ titj cos {Oi -Oj)- titj + A ^ ti , (34) 

{ij) {ij) i 

where the first two sums run over nearest-neighbor pairs and the last one is over all lattice 
sites, except those at the surface. In this lattice gas model of ^He-^He binary mixtures the 
coupling constant K and the field A are related to the effective "'He-^He interactions Kap 



(see, e.g., Ref. [23]), 

K = K33 + Ku - (35) 
and to the chemical potentials fi^ and /i4 of ^He and ^He, respectively, 

A = ^3 - /i4 + 2g(i^33 - ^^34) , (36) 

where q is the coordination number of the lattice {q = 2d, where d is the spatial dimension 
of the system; g = 6 in the present case). In the liquid the effective interactions Ka/3 are 
different for different a and (3 due to the differences in mass and of statistics between ^He 
and ^He atoms. 
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The properties of the model described by the bulk Hamiltonian Hb have been studied 
within MFT and by Monte Carlo simulations in = 3 [20] . In contrast to its two-dimensional 
version, for which there is no true tricritical point for any value of the model parameters, 
in (i = 3 for reasonable values of the interaction parameters the resulting phase diagram 
resembles that observed experimentally for ^He-'^He mixtures, for which phase separation 
occurs as a consequence of the superfluid transition (see Fig. H]). The form of the surface 
Hamiltonian Tig should capture the phenomenon of superfluid film formation near a wall in 
^He-'^He mixtures 18| which generates an effective repulsion of ^He atoms by the wall. The 
van der Waals interactions between the wall and ^He or ^He atoms are equal. However, ^He 
atoms occupy a larger volume because of their larger zero-point motions. This gives rise to 
the preferential adsorption of ^He atoms at the substrate-fiuid interface, which may induce 
a local superfluid ordering and an enrichment of ^He near the opposing fluid-vapor interface. 
Here we choose the following form for Tig: 

K = 5A« J]t, + <5AM^t,, (37) 

i i 

where the first sum runs over the sites of the first layer and the second over those in the L-th 
layer of the lattice. The differences SA^^^ = A*^'^ — A and 6A^^^ = A^^^ — A are measures of 
the relative preferences of "^He atoms for the two surfaces such that (5A^') < corresponds 
to the preference of "^He atoms for the solid substrate. 

B. Mean-Field Theory 

We have studied the above model for the film geometry within mean-field theory. We 
have employed the variational method based on approximating the total equilibrium density 



distribution by a product of local site densities pi (see, e.g., Ref. [28|). The corresponding 
variation theorem for the free energy reads 

F<F, = TiipH) + (l//3)Tr(plnp), (38) 

where F is the exact free energy and Fp is an approximate free energy associated with the 
density distribution p; P = l/^ksT). The minimum of Fp with respect to p subject to the 
constraint Trp = 1 is attained for the equilibrium density distribution p = e~^^/Tr(e^^^). 
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Within mean-field theory the density distribution in the film geometry is approximated by 

L 

p^Po^A'^Hpi, (39) 

i=l 

i.e., the density distribution is constant within each layer parallel to the surfaces but varies 
from layer to layer. We treat the local layer density pi as a variational ansatz, and the best 
functional form in terms of ti and 9i is obtained by minimizing Fp^^/ A + r]Tr{pi) with respect 
to Pi and with 77 as a Lagrange multiplier in order to implement Trp = 1. This leads to 

Pi = e-^'*VTr(e-^'^0 (40) 
where hi is the single-layer mean field given by 

hi= - J{Mt\ + q\\M^'^ + M^l\)t, cose, -J {Ml% + q\\M^^^ +mII^^^^^ 

- K{Qi_, + q\\Qi + Qi+^)ti + A«t,, (41) 

where A^*) = A for i 7^ 1,Z, and A'^*) = A*^^)(A'^'')) for i = 1{L). We have introduced the 
following order parameters: 

Qi = l- X{i) = Tr{tiPi) (42) 

and 

M^'^ = Tr {piti cos Oi), Mf ^ = Tr (p^t, sin Oi). (43) 

Qi — {ti) corresponds to the concentration profile of ^He, X{i) — l — {ti) to the concentration 
profile of ^He, and m/^^ , M-^^^ are the components of the two-component superfiuid OP profile 
Mj . g| I is the in-layer coordination number while each site (but not in the first and last layer) 
is connected to q' atoms in each adjacent layer and q = qw + 'iq' is the coordination number in 
the bulk of the lattice. Within our model q' = 1 and q\\ = 2{d— 1). This yields the following 
set of self-consistent equations for the OP Mj = {M^^\ 0) = (m^, 0) in the ith layer : 

Qi = WJbi)/ (^e-^(^"-^''^) + WJk)) , (44) 

and 

= /i(/5J6.)/ (e-'^^^"^-^'*^) + /o(/3J6.)) • (45) 

Io{z) and h{z) are the modified Bessel functions of the first kind, T is the temperature. We 
have introduced 

bi = rui-i + q\\mi + rrii+i for i y^l,L, (46) 
15 



bi = q\\mi + 1712, and bi = mi_^i + q\\mi, and analogously 

ai = Qi„i + q\\Qi + Qi+i for i y^l,L, (47) 

o-i = + Q2, and flj^ = + qwQi- The coupled sets of equations for Qi and are 
solved numerically by standard methods of multidimensional root finding. The equilibrium 
solution minimizes the free energy per number of lateral lattice sites JF = Fp^JA: 

L-l r 



i=2 



J K 

-{mi-ivfii + q\\m1 + mi+inii) + —{Qi-iQi + qwQl + Qi+iQi) 



+ fcBT^ln(l-Q,) + /i + /2, (48) 



where 

J K 

fi = 2 (^11"^? + "^2rni) + y (gilQ? + (49) 

and 

J K 

/2 = -^{.rni_^mi + g||m|) + —{Qi_iQi + (50) 

The above equations neglect the helicity, i.e., Mj = {M^^\ 0) = (mj, 0). In general the helic- 
ity might be non-zero because the BC for the superfiuid OP are effectively non-symmetric, 
i.e.. Ml 7^ whereas M/, = so that the superfiuid OP can in principle rotate across the 
film. The relevance of the helicity on the Casimir force will be analyzed elsewhere. 



C. Results for ^He-^He mixtures 

First, we have analyzed the semi- infinite system. Close to the line of hulk critical points 
we have found a higher ^He concentration near the surface (chosen to be the left side of the 
system), which induces a local superfiuid ordering. By varying T and A one obtains a line of 
continuous surface transitions corresponding to the onset of the formation of this superfiuid 
film near the wall; it meets the line of bulk critical points at a so-called special transition 
point, the position of which depends on the value of A*^') (see Fig. H]). These findings are in 



agreement with the results of a Migdal-Kadanoff analysis [29 1. 

In the film geometry the Casimir force fc (Eq. ([2])) is obtained by calculating f^^{L) 
(see Eq. (jS])) for L and L + 1 and taking the difference. (Note that in the lattice model / 
is the total free energy of the film per number LA of lattice sites and fj, is the bulk free 
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energy density per LA. Accordingly f^^iL) = (/ — fb)L/ (kBTt), fc = —df^^/dL, as well as 
i) = fcL'^ with d = 3 near tricriticality and d = 4 near the A-transition are dimensionless. 
In order to avoid a clumsy notation we do not introduce different symbols for the lattice and 
the continuum versions of the free energies.) Figure summarizes our result for a film of 
thickness L = 20, K/J = 0.5, A^V-^^ = -3, and A^^) = At/ J ~ 0.61, which is the tricritical 
bulk value. Such a choice of the surface coupling constants corresponds to non-symmetric BC 
and is consistent with the assumption made in Ref. Q| for the concentration profile across the 
wetting film, whereupon at the interface with the vapor the ^He concentration takes the bulk 
value. For temperatures above the bulk coexistence line at first-order demixing transitions 
fc is calculated along the thermodynamic paths indicated in Fig. H] which correspond to 
fixed '^He concentrations X. Our selection of X covers the tricritical region as well as the 
crossover to the critical superfiuid behavior of pure "^He, i.e., X = 0. In order to calculate 
the force at a fixed value Xq we first determine A(X = Xq,T) by solving the two coupled 
self-consistent equations for Q{A,T) = 1 — X and M(A,T) in the bulk (Eqs. (12) and (13) 
in Ref. 20[). For each temperature along the thermodynamic paths indicated in Fig. Hlwe 
solve Eqs. f H4l) and ( H5l) with this value A(X = Xq,T). This renders the profiles Q{1) and 
m{l) and allows us to calculate the free energy from Eq. fHHj) . When upon lowering the 
temperature the paths of constant X reach the coexistence line of two-phase coexistence 
(see Fig. Hj) we continue our calculations along the coexistence line, infinitesimally on the 
superfiuid branch of bulk coexistence. In Fig. [5] this leads to the full line for T < Tf, i.e., 
y<0. 

Contrary to the LG model, for the present microscopic model it is natural to express 
the properties of the system as functions of the experimental thermodynamic fields t and 
(A — At)/{kBTt) or the scaling fields t and g (see Eq. flTTl)). Accordingly, we present our 
results for the Casimir force in terms of the scaling function defined through the relation 
■(9 = L^fc as a function of only a single scaling variable y = tL^^'^ = {{L/ a) / {^l^^))^/" . 
^0 = is the amplitude of the order parameter correlation length ^ = ^^t^'^ = above 
Tt and z/((i = 3) = 1. The second relevant scaling variable x = gL? also varies along a 
path of fixed ^He concentration (see Fig. [6]) and a proper scaling description has to account 
for it. However, in order to be able to compare our results with the presentation of the 
corresponding experimental ones [4], we follow Ref. |^ where the variation of x has been 
neglected. As can be inferred from the phase diagram in Fig. El the (yf-components of the 
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paths X = const in the phase diagram are smaller than the t-components, so that the form of 
the scaling function for these paths are expected be close to 'd{x = 0,y). Also experimentally 
the variation of the scaling variable g along the path of fixed X cannot be determined easily. 

Near the tricritical point paths of constant X cross three different phase transition lines: 
the surface transition line, the line of bulk critical points, and the line of first-order phase 
coexistence. As shown in Fig. [5], close to the surface transition fc is small and this transition 
does not leave a visible trace in its behavior, fc remains small up to the coexistence line 
or to the line of bulk critical points for X > Xf or X < X^, respectively. There it increases 
very steeply and for ^He concentrations X < Xt upon crossing the line of bulk critical points 
there is a break in slope (see the dots in Fig. giving rise to the formation of shoulders 
which are similar to those observed experimentally |4j]. When T reaches the temperature of 
first-order phase separation, fc is given by the curve (full line for ?/ < in Fig. [5]) common 
to all values of X. These curves of constant X meet the full line with different slopes. 

The aforementioned common curve exhibits a pronounced maximum below at ?/ ~ 
—0.74 and gradually decreases to zero for y — oo. The properties of the Casimir force in 
this temperature region can be attributed to purely interfacial effects. Indeed, we observe 
that below Tf both the concentration and the superfiuid OP profile corresponding to this 
common curve display an interface-like structure separating two domains of the coexisting 
bulk phases (see the case t = —0.0633 in Fig. [7]). This film phase is soft with respect 
to shifts of the interface position and is similar to the one occurring in Ising-like films 



with opposite BC 



30| for temperatures below the bulk critical temperature but above the 



wetting temperature of the confining walls, in which case the Casimir force is repulsive 
with a pronounced maximum occurring 

a positive sign of the force can be regarded as a consequence of entropic repulsion 32|. 
Typically the maximum of the force occurs at that temperature T at which the interfacial 
width, which is proportional to the bulk correlation length ^ of the order parameter, becomes 
comparable with the width L of the film. In the present case both the concentration and the 
superfiuid OP profile contribute to the free energy and hence to the Casimir force. Their 
interfacial widths are proportional to correlation length ( associated with concentration 
fluctuations and to the OP correlation length ^, respectively. As can be seen from Figs.[7]and 
[HI within MFT these interfacial widths and therefore ( and C, are comparable. Accordingly, 
by analogy with Ising-like systems 301] we expect that within MFT the maximum of the 
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force occurs when ^ (or, equivalently, C — is of the order of L, which is actually consistent 
with what is observed in Fig. [5], where the maximum of the scaling function is located at 
?/ ~ — 1. We may expect that also in the actual system the occurrence of the maximum 
of the Casimir force below the tricritical point can be attributed to such interfacial effects. 
However, since the correlation length of the superfluid OP ^ = oo in the superfluid phase it 
is not yet clear which length scale governs the interfacial width of the superfluid OP profile 
in the 'soft mode' phase below Tt and hence what length scale determines the position of 
the force maximum. 

For X < Xt — 0.05 we observe a crossover to the critical superfluid behavior of pure 
^He and a gradual formation of a second, less pronounced local maximum located slightly 
below the line of bulk critical points (?/ > in Fig. This local maximum decreases upon 
departing from Xt and finally fc becomes vanishingly small along paths which cross the 
line of bulk critical points above the special transition S (see Fig. Hj). This is expected, 
because above S there is no longer a superfluid film formation near the solid substrate for 
thermodynamic states corresponding to the bulk "normal" phase of a fluid close to the 
line of bulk critical points. This means that the superfluid OP in the film is identically 
zero up to the line of bulk critical points and the BC effectively turn into the type (O, O) 
for which fc vanishes within MFT. (For (0,0) BC fluctuations beyond MFT generate an 
attractive Casimir force /c < [l^.) For lower T, fc increases steeply upon approaching 
bulk coexistence revealing that interfacial effects associated with the 'soft mode' lead to 
a much stronger Casimir effect than the critical fluctuations near the line of bulk critical 
points. 

IV. RESULTS FOR PURE ^HE 

A. The limiting case of the VBEG model 

In this section we consider the limiting case A — > —oo in which all lattice sites are 
occupied, i.e., ti 1. In this case the first term of the bulk Hamiltonian Tib in Eq- dM]) 
corresponds to the classical XY model (the planar rotator model) for pure ^He and therefore, 
as far as the bulk contribution is concerned, the partition function of the VBEG model 
reduces to that of the XY model up to a factor e^^^ where N is the number of lattice sites. 
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The corresponding MFT equations for the bulk OP can be inferred from Eqs. (jHj) and (H5|) 
with rrii = M yielding 

for temperatures below the bulk superfluid transition, which is located at Ts{X = 0) = 
T\ = qJ/2, and Q = 1, M = above T\ = Ts{X = 0). The scaling behavior of the free 
energy and of the Casimir force close to this critical point (see below) is consistent with 
an upper critical spatial dimension d* = 4. The crossover to the tricritical behavior with 
d* = 3 and with tricritical exponents occurs only upon approaching the tricritical point 
A = (Tt/T,(0) = 2/3, Xt = 1/3) (see Fig.SD. 

In the slab geometry we take also the limits A*^'-*, A*^''-' — — oo which, together with the 
absence of external fields coupling to the superfluid OP, lead to (O, O) BC for the superfluid 
OP. Thus this limiting case allows us to study the Casimir force for wetting films of pure 
'^He near the superfluid transition at Tc = Tx. We remark that in the slab geometry the 
superfluid transition is actually of the Kosterlitz-Thouless type 33|| . However, this change of 
the character of the transition is not captured by MFT. The corresponding set of equations 
for the superfluid OP in the Z-th layer of the slab is: 

m = 'i^,, k^m,_, + g,m, + m,^, for l^l,L, (52) 

where bi = q\\mi + m2 and bi = mi_i + q\\mi. The equilibrium free energy divided by the 
number A of lattice sites within one layer takes the form 



L-l 



+ + m2mi) + ^{mi_imi + q\\ml) 



1=2 

+ knTj^HWJbi)) (53) 

1=1 

where m/, / = 1, . . . , L, are the solutions of Eq. flS21) . Solving Eq. (1521) for different widths of 
the film we have found that the superfluid OP profile vanishes for temperatures larger than 
a certain Tc{L) < Ts{X = 0) = Tx which can be identified with the critical temperature 
Tc{L) of the slabs. Below Tc{L) the corresponding Casimir force turns out to be negative 
(i.e., attractive) as expected for (0,0) BC pertinent to the case of pure "^He. The lattice 
calculations have been carried out for = 3 and are presented in terms of the scaling function 
^^{y = r(L/^+)2) = L'^fc with = 4 in accordance with MFT and r = {T-Tx)/Tx. Within 
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lattice MFT the actual space dimensionality d of the lattice does not influence the shape of 
the scaling function in the scaling limit ^(r > 0) = ^qT^" ^ a; indeed, it enters only into the 
non-universal amplitude C,q via the ratio q'/q = (2d)~^ between the bulk inter-layer and the 
total site coordination numbers q' and q, respectively, has been calculated for L = 20, 40, 
and 60 and is plotted in Fig. M function of y = T{L/^^y^^ with the MFT value z/ = 1/2. 
Exploiting the fact that within MFT ^(r < 0) is flnite we have determined the amplitude 
of the correlation length ^{t < 0) = ^o{—t)^'^ from the exponential approach of the OP 
profiles towards the corresponding bulk values ruf, which are actually attained in the middle 
of the film (see Fig. [10]) at temperatures sufficiently below Ts{X = 0) (see, e.g.. Fig. 21 



in Ref. 



34(1). The MFT universal amplitude ratio /^q = V2 then yields the estimate 



C,^ — 0.41 for the VBEG model on the lattice. We emphasize here that scaling of the force 
data occurs only for surprisingly thick films, i.e., L > 60, as revealed clearly by the analysis 
presented in the next subsection. 



B. Comparison with the Landau- Ginzburg theory 



In Ref. 



34| within MFT for the 0(2) LG continuum theory (see Eq. ([T]) with vq = 0) 



the order parameter profiles $ = {m{z),0) in a slab with (O, O) BC have been calculated 
analytically (see Eqs. (202) and (203) in Appendix D in Ref. [34)]). It turns out that as 
a function of the scaling variable y = t{L/E,q Y^'^ = tq-L^ (where tq oc r is the coefficient 
appearing in Eq. ([1])) the mean-field OP profile m{z) vanishes for y > ym = ~7r^, whereas it 
is nontrivial for y < y^, breaking the original 0(2) symmetry. This occurs for temperatures 



below t 



le shifted critical point of the film which therefore corresponds to y = ym (see 



Ref. I35]). In Fig. [10] we compare the OP profiles (normalized by the corresponding bulk 
values as to obtain universal scaling functions of y and z/L) calculated within the VBEG 
model (for a lattice with L = 150) and within LG continuum theory for a selection of values 
of the scaling variable y. The agreement between the profiles is very good, although the 
VBEG profiles exhibit a slight asymmetry with respect to z/L = 1/2 which is due to the 
limited numerical accuracy of the lattice calculation. 

The knowledge of the analytic expression for m{z) allows one to compute the stress tensor 
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(Eq. (jlj)) as a function of the scaling variable y: 

T„ = = 0))^ = UJlTW Uj ■ ^ < ^" = ■ (54) 

[ 0, for y>ym, 

where Am = 37r^/(2'Uo) and k = k{y < ym) is the real solution of the implicit equation 

JL = ±(i + k^)K\k) (55) 

ym TT^ 

where K{k) is the complete elliptic integral of the first kind, such that k{y = ym) = and 
k{y — > — oo) = 1. The stress tensor %z,b in the bulk, related to the bulk free energy density 
/^(r), can be obtained from Eq. flM|) in the limit L ^ oo at fixed reduced temperature r, 
yielding T^z^bij < 0) = AmL'^^y/ymY (which is actually independent of L due to ?/ oc rL^) 
and Tzz,b{T > 0) = 0. Accordingly, the Casimir force fc per unit area of the cross section of 
the film and in units of ksTx is given by fc = Tzz — Tzz,b and its scaling function 'Oq^ = L'^fc 
can be derived from the expressions for Tzz and Tzz,b discussed above: 

2\ 2 / „, \ 2 



<^(l/) 



1-PV / y 



im 

' 2 



-Am\^] for ym<y<^ ^^^^ 
for ?/ > 0. 



The independent calculation of ^^^{y)-, recently presented in Ref. 36|, agrees with this 
expression. At y = ym = — tt^ the scaling function (!56l) exhibits a cusp singularity at which 
it attains its minimum value t?^^ = ^o'~^{ym) = —Am < where Am is given after Eq. fl5^ . 
Within MET the coupling constant Uq and therefore Am remain undetermined. In order to 
compare the LG result with the VBEG results, accounting also for corrections due to the 
finite size L of the latter, we introduce an adjusted scaling function ^Q^{y) which is given 
by Eq. (156!) with Am = Am{L) and ym = ym{L) determined by a best fit to the VBEG 
scaling function '&o{y) calculated for lattices with L = 20, 40, and 60. For all values of L 
considered, ^Q^{y) provides a very good fit to the numerical data, as demonstrated in Fig. [H] 
for L = 60. In the inset of Fig. [H] we plot the functions Am{L) and ym{L) obtained from 
the fit. According to the results of the LG theory one expects ym{L — > oo) = — tt^ ~ —9.87 
(which is represented as a solid line in the inset), and indeed the results of the VBEG model 
show the correct trend, although finite-L corrections are still present even for the largest 
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lattice L = 60 considered here, with ym{L = 60) ^ —9.31. The amphtude Ajn{L) shows 
even stronger corrections and indeed the value Am{L = 60) ~ 2.45 might underestimate the 
actual asymptotic value by 15-20%. 

Beyond MFT the renormalized coupling constant u attains its fixed-point value under 
RG fiow which fixes the amplitude Am and the magnitude of the corrections to the scaling 
functions. This would then allow a complete numerical test with the scaling function -t^o of 



the VBEG model as obtained, e.g., from Monte Carlo simulations. In Ref. [36| the amplitude 
Am = 37r^/(2Mo) (see the text after Eq. fl5^ ) has been estimated beyond MFT by replacing 
Mo by the fixed-point value u* calculated within field theory. Although this approach provides 
a theoretical estimate for {Am)theo = 6.92, it fails in accounting quantitatively for the actual 
amplitude {Am)exp = 1-30 ± 0.03 observed in experiments [3]. 

For a given film thickness L, the position of the minimum of the scaling function cor- 
responds to the reduced critical temperature Tm{L) = {Tc{L) — Tx)/Tx = Umiit / LY^" 
which reflects the onset temperature Tc(L) < Tc(L = oo) = Tx for superfluidity in the 
slab. For t > Tm the superfluid OP proflle vanishes and so does the mean-fleld free en- 
ergy of the fllm. Thus from Eqs. ^ and ([3]) it follows that for T > Tc{L) one has 
L'^/c = -L'^fh/ikBTx) ~ -L'^r^-" = -{tL^/^Y", using the hyperscaling relation 2-a = du. 
For d = A and u = 1/2 this implies that 'do{ym < y < 0) ^ y'^ (for y > 0, within MFT fb = 
and therefore fc = 0) which agrees with Eq. fl56l) . 

V. DISCUSSION OF THE RESULTS OBTAINED FROM THE VBEG MODEL 
A. '^He-^He mixtures 

As one can infer from the comparison of Figs. E] and [3] the qualitative features of the 
scaling functions for ^He-^He mixtures extracted from the experimental data for X ^ Xt, 
such as the sign of the force, the occurrence of the pronounced maximum below Tj, and the 
formation of shoulders above Tj, are well captured by the present lattice model. The breaks 
in slopes upon crossing the A-line shown in Fig. [5] are features of the mean-fleld approach 
and expected to be smeared out by fluctuations. 

The experimental data for the Casimir force fc exhibit a maximum at tL ~ — IsA which 
cannot be related to the condition ^ ~ ~ L borne out by the mean-fleld analysis (Fig. [5]) 
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because actually, i.e., beyond MFT, ^ = oo in the superfluid phase. Further studies are 
needed to determine what length scale governs the interfacial width of the superfluid OP 
profile in the 'soft mode' phase below Tf. This analysis, which is left to future research, has 
to take into account that the actual width of the interface formed in the film (see, e.g., the 
case if: = —0.0633 in Fig. [7]), is broadened both by the Goldstone modes in the superfluid 
phase and by capillary-wave like fluctuation. 

Different from the mean-field scaling function the experimental one does not vanish at 
low temperatures, which is expected to be due to the aforementioned Goldstone modes of the 
broken continuous symmetry in the superfluid phase and due to helium-specific [l5| surface 
induced fluctuations which both evade the present mean-field analysis. A similar behavior 
has been found in wetting experiments for pure ^He films near the A-line j2|, in which the 
film thicknesses above and below the A transition are not the same, so that the wetting films 



are thinner in the superfluid phase. For pure ^He Zandi et al. 15| pointed out that the 
Goldstone modes indeed lead to thinner superfluid films for T ^ Tc. But this estimate is 
not applicable for T T\ and for T ^ Ta it is too small to account for the experimentally 
observed magnitude of the thinning. This view of the effect of the Goldstone modes on d 
is supported by Monte Carlo simulation data for the XY model with periodic BC [3^. The 
capillary wavelike surface fluctuations, which occur on one of the bounding surface of the 
superfluid ^He wetting film, give rise to an additional force (similar in form but larger in 
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magnitude) which may then together explain the experimental observation |l5l . 

For a mixture, however, it is possible that the apparent thickening of a wetting film 
as inferred from capacity measurements might be, at least partially, an artifact due to a 
significant change of the permittivity within the film [39.] . Upon inferring the film thickness 
from the permittivity, in Ref. it was assumed that Xjum = Xt which does not hold 
at low temperatures at which the 'soft mode' occurs. In order to estimate the error the 
assumption XfUm = Xt introduces into the determination of the film thickness L we repeat 
the calculation for determining L by taking into account the interface-like concentration 
profile below Tt (see Fig. [7]) and by assuming a mean field-like shape: 

X{z) = 1{X' + X") - ^(X^ - tanh[(z - zo)/(2C)], (57) 

where X^ and X^^ are the concentrations of the coexisting bulk phases (see the triangle 
in Fig. H]), zq = L/2 is the position of the center of the interface, and ( is the correlation 
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length associated with concentration fluctuations. We note that ( is finite in the superfiuid 
phase whereas .^ = 00 for the superfiuid OP. The effective permittivity constant efum of the 
film follows from adding in series the capacitance C for each slice of the film and from using 



^faUX,T)= (58) 
Jo dz/fiz) 



where e{z) is related to the concentration profile via [40] 



e{z) - 1 = (5.697 - 1.402X(^)) x 10"^ (59) 

From this we have found that neglecting at low temperatures the variation of the con- 
centration across the film introduces an error in the determination of its thickness from 
capacity measurements (leading indeed to an increased film thickness) which is about 35% 
of the 40A difference in thickness reported above and below Tt. Specifically, at T = 0.65K 
the bulk concentrations are = 0.325, X^^ = 0.825, and the bulk correlation length 
is C = Col^l ""^ ~ 5.1A, where following Ref. [4] we have assumed Co = 1.3A as the value 
measured for concentration fiuctuations far below Tt in the superfiuid phase. Accordingly, 
approximating the actual inhomogeneous permittivity by the homogeneous one gives rise to 
an error ^ I4A. 

In the crossover regime along the line of critical points connecting the tricritical point and 
the critical A-transition in pure "^He only few experimental data for the thicknesses of the 
wetting films are published. Nonetheless, the observed variations of film thicknesses there 
again agree with the present theoretical findings for the Casimir force. In particular, one 
observes a rapid thickening of the films upon approaching the line of bulk critical points; for 
specific values of X a small maximum located slightly below the line of bulk critical points 
is also visible (compare Fig. E]). 

Two reasons impede a more quantitative comparison of our results obtained within the 
VBEG model with the experimental ones. First, for our choice of surface terms in the 
Hamiltonian the fixed-point BC (+, O) for the order parameter cannot be realized within 
the VBEG model. Taking the limits Ai — >■ — cx) and A2 — 00 in Eqs. f l44p and 0451) assures 
that X{1) = and X{L) = 1. However, even this limiting concentration profile does not 
induce the required BC: although m(L) = one has m(l) = /i(/9J6i)//o(/3J&i) 7^ 1, i.e., 
the superfiuid OP at the solid substrate is never saturated at its maximum value 1 which 
corresponds to the BC (-I-) (see Fig. [8]). We have checked that in this limiting case with 
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respect to Ai 2 the qualitative behavior of the Casimir force is the same; only the magnitude 
of fc is slightly bigger ("(9(0) ~ 0.5 for the limiting case, whereas 19(0) ~ 0.4 for the case 
shown in Fig. In order to be able to extract universal properties - which requires to 
reach the fixed-point BC - it would be necessary to introduce a surface field which couples 
directly to the superfiuid OP so that the BC (+) can be realized; but such a surface field 
has no physical basis. Finally, even at the upper critical dimension d = d* = 3 due to 
logarithmic corrections our present MFT is not sufficient. However, a naive correction of 
obtained within the VBEG model by multiplying it by the logarithmic factor (ln(L//o))^^^ 
(see Eq. ( 133|) ) derived within the LG model does not capture the proper universal scaling 
behavior. Instead renormalization group schemes for the VBEG model have to be employed. 

Nonetheless, our MFT results for the scaling function ?9 within the VBEG model and 
for X = Xt, a matched with respect to its amplitude with the experimental data at the 
tricritical point y = and after adjusting the scaling variable y hj a factor yth such that 
the experimental and theoretical positions of the maximum of the scaling function are the 
same (which is achieved for yth ~ 0.065), yield an adjusted scaling function i){y) which 
reproduces rather well the experimental data (see Fig. [TTi) . especially near the maximum 
where interfacial effects are expected to be dominant. This observation is consistent with 
our interpretation that the formation of this maximum is dominated by the occurrence of 
the 'soft mode' phase which does not depend on the details of the surface fields. We note 
that according to Fig. [11] the experimental data nominally for X = Xt more closely agree 
with the theoretical ones for X = Xt — 0.01. This raises the question as to whether the 
experimental ^He concentration in the film is actually shifted relative to the bulk one. 

B. Pure ^He 

The theoretical models discussed in the previous sections (VBEG and LG as lattice and 
continuum models, respectively) capture the universal features of the collective behavior 
close to critical (and tricritical) points, such as the Casimir force. (These models have 
no predictive power concerning non-universal properties.) The associated finite-size scaling 
functions acquire universal forms if expressed in terms of proper scaling variables, such as 
L/^(r), where ^(r) is the correlation length which controls the large-distance exponential 
decay of the two-point correlation functions of the OP fluctuations in the bulk at the reduced 
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temperature r. In systems with discrete symmetry one has ^(t ^ ) = Co ("''") ^ 



^,{r ^ 0+) = Co r^'', where Co are non-universal, i.e., system-dependent, amphtudes such 
that the ratio Co^/C(7 is universal (see, e.g., Ref. 4l|)- Accordingly, the scaling function 
maintains its universal character also as a function of y = t^L/^qY^'^ in the notation of 
Sec. llVl or. alternatively, t(L/^qY^'^ . However, in the case of pure ^He, the bulk correlation 
length C(t < 0) below the A-transition is infinite due to Goldstone modes and therefore 
C(7 cannot be defined directly from the behavior of C(t < 0). Alternatively, one might 
define a different length scale ^'^{r < 0) = C,o{—t)^'^ associated with the power-law decay of 
tranverse correlations in the superfiuid phase, which is related to the superfiuid density; the 



non- universal amplitude C(f forms a universal ratio with Cfj" (see, e.g., Refs. 4l|, |42|). For pure 



'^He, experimental estimates of (C(f)ea;p range from 1.2A 43] to 3.6A 44], depending on the 
way it is measured. In view of this experimental uncertainty and of the complication related 
to the introduction of C(f oc ^/uo^Q [42] within the MFT discussed in Sec. HVt we present 
the comparison between experimental data and the VBEG model in terms of the scaling 
variable y, which involves the non-universal amplitude Co^ the value of which is well assessed 
experimentally for ^He, (Co')ea,'p = 1.43A at saturated vapor pressure 45[, and theoretically 
for the VBEG model, Co" = 0.41a within the present MFT, where a is the lattice spacing 
(see the end of Subsec. IIV A|) . Within the LG model one has an analytic expression for 
in terms of the parameters of the model (see Eq. (6.4) in Ref. 13| for obtained within 
the dimensional regularization scheme). 

In Fig. [12] we compare the scaling function obtained from the experimental data for the 
case of pure ^He ^ (for a film thickness L = 423 A 39]) with the MF scaling function i)o{y) 
of the VBEG model which is universal for sufficiently thick films. The scaling functions 
are normalized by their absolute values ]'(9mm] at the minimum. In order to summarize all 



available theoretical results we report in the right inset of Fig. [12] the comparison between 
the experimental data for T > T\ and the scaling function obtained from the field-theoretical 
e-expansion (e = 4 — ci) [13]] as follows: The scaling function 6+0,0(2/+) of the finite-size 
contribution of the renormalized free energy / provided in Eq. (6.12) of Ref. 3] has been re- 
expressed for = 2 (XY model) as a function of ?/ = T{L/^'^y/'' viay^ = y^^'^{l+e/lQhiy) + 
O(e^) (where ?/+ is defined after Eq. (4.6) in Ref. [3]). The resulting expression 6+0,0(2/) = 
^0(2/) + ^0i{y) + 0{e^) is then extrapolated to three dimensions e = 1 either as Q^looin) = 
6o{y) + e^i(2/) (yielding the solid line in the inset) or <d%]o{y) = 0o{y)/[^ - ^0i{y)/6o{y)] 
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(dashed line), corresponding to the Pade approximants [1,0] and [0,1]. The scahng function 
of the Casimir force is then provided by 'diy) = {d — l)6+,o,o(2/) ~ iv / ^)^'+ o oiv) where 
d = 3 and z/ ~ 0.67 (see, e.g.. Table 19 in Ref. |4l)]), accounting for the actual expression of 
the scaling variable y in three dimensions. 

Discrepancies, such as the position of the minimum, the shape of the scaling func- 
tion for y > ym, the behavior for y — > — oo, and the nonvanishing of -^^^^ for y > can 
be attributed to fluctuation effects neglected in the present MF approach. Field-theoretic 
renormalization group calculations beyond MFT yield a quantitative agreement with the 



experimental data for y > 10, |l3|, |lj] (see Fig. [T2l) : however, so far this field-theoretical 



approach cannot be extended to the case y < [38]. From the analysis of Subsec. IIVBI 
it follows that for fixed L the position y^ = — vr^ ~ —9.87 of the minimum is associated 
with the critical temperature Tc{L) of the film. The experimental data exhibit the posi- 



tion of the minimum at Xr 



9.8 ± 0.8A^^^ where x = tL^^' 



corresponding 



to {ym)exp = Xmin/{^o)exp — — 5.7 ± 0.5 which is consistent with the experimental indi- 
cation in the sense that the onset of superfluidity in the films occurs within the range 
— 12A^^'' ^ X < —ih^^^ [7], i.e., — 8 < ?/ < —5. But these values of y are considerably 
larger than the value — vr^ predicted by the LG approximation. In spite of the shortcomings 
mentioned above the comparison between the experimental and theoretical scaling func- 
tion is nonetheless encouraging. The present MF approach does not address the issue that 



\'drmnim\ 



exp 



20 0, 0] whereas theoretically this ratio is ~ 1 for periodic EC j38| ; it is 
difficult to expect that this ratio reaches the experimental value 20 corresponding to the 
actual (O, O) BC. 

In passing we mention that in Ref. 36| the comparison between Eq. ( !56|) and the ex- 
perimental data of Refs. 2|, |7|] is seemingly affected by an inconsistent normalization of the 
experimental and theoretical scaling functions which are actually plotted as a function of 
T{L/^^y/'' (with taken from Ref. |43]) and t(L/^qY^'^ , respectively. This artificially 
reduces the resulting discrepancy between the experimental and theoretical results in com- 
parison to the one displayed in Fig. [T2i 
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VI. SUMMARY AND OUTLOOK 



Based on mean-field analyses of the vectoralized Blume-Emery-Griffiths model and of the 
continuum Landau- Ginzburg theory as well as by applying renormalization group analyses 
we have obtained the following main results: 

(1) By using mean- field theory, near the tri critical point (Fig. Hj) we have calculated 
the scaling functions of the Casimir force within the continuum Landau-Ginzburg theory 
(Eq. (Q) for the 0(2) model of ^He-^He films of thickness L (see Figs. [T]and[2]). The scaling 
functions depend on two relevant scaling variables Uq and vq (see Eq. f|T8|) ). By fitting the 
amplitude of the scaling variable and the amplitude of the Casimir force, which remains 
undetermined within the LG mean-field approach, one finds a reasonable agreement with 
the experimental data along the thermodynamic path of constant tricritical concentration 
of ^He (see Fig. ED . 

(2) The application of fieldtheoretic renormalization group analysis in spatial dimension 
d = 3 yields the correct asymptotic leading behavior of the Casimir force at the tricritical 
point. As a function of the film thickness L it has the form of a power law ~ L^^ multiplied 
by the square root of the logarithm of L and by the universal Casimir amplitude (Eq. 0231) ). 

(3) Using the fieldtheoretic renormalization group analysis we have derived the form of 
the finite-size scaling for the Casimir force in the vicinity of the tricritical point and have 
obtained renormalized mean field scaling functions (see Figs. [T]and[2]). It turns out that 
also one of the scaling variables acquires a logarithmic correction (Eq. ( l33i) ). 

(4) Using mean-field approximation we have calculated the scaling function of the Casimir 
force within the vectoralized Blume-Emery-Griffith lattice model of ^He-^He mixtures along 
the thermodynamic paths of fixed ^He concentrations (see Figs. HI O and [5]). For concentra- 
tions of ^He close to the tricritical concentration our results are in a qualitative agreement 
with the available experimental data (see Figs. [3|, O and [11]). Our calculations also predict 
the crossover behavior of the Casimir force along the line of critical points connecting the 
tricritical point and the A-transition for pure ^He. We have found that the pronounced 
maximum of the Casimir force, which occurs below the tricritical temperature, is associated 
with the formation of a 'soft mode' phase within the film (see Figs. [7|and[8]). 

(5) We have analyzed the limiting case of the VBEG model which corresponds to the 
classical XY model for pure "^He. Within mean-field theory we have been able to show 
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that for sufficiently thick ffims the scahng functions as obtained from the lattice model for 
the Casimir force are in an agreement with the ones obtained from the continuum 0(2) 
Landau-Ginzburg theory (see Fig. [H]). The encouraging comparison of the former with the 
experimental data is displayed in Fig. [T2i 

As an outlook we propose to test experimentally the scaling of the Casimir force for 
different thicknesses of the wetting films by taking into account logarithmic corrections. 
Moreover it appears to be promising to study experimentally in more detail the crossover 
of the Casimir forces between their tricritical behavior and their critical behavior near the 
A-transition and to compare it with the theoretical predictions presented here. 

A.M. benefited from discussions with R. Garcia and M. Krech. 
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FIG. 1: Dimensionless MF scaling function ^^^{roL'^ = 0,y*'^^) = /ci^(wo/90)^/2 (ggg gq. (fTo]l ) 
with y^"'^^ = (5/(2vo))^^^'Uo-^ ~ tL plotted together with the renormalized mean field scaling func- 
tion fcL^ = 7?^*^^(0,y^^-^^) (see Eq. 1^ and the main text) with y^^^^ = uL{\n{L/lo)f/^'^ , 
u = 7u/(247r2), and L/Iq = 400. T?*^^(0,y*^^ ^ cx)) ~ (thin dash-dotted line) 

and ^^^^{0,y^'^^ — > 0) ~ 2.76 — 0.605y*^^ (thin dashed line). The asymptotic behavior of 
^RMF ^Q^yRMF^ can be obtained from the one of -d^^^ {0,y^^) by multiplying the ordinate by 
the factor (28/(87rV3))^/^(ln(L//o))^/^ and the abscissa by the factor (ln(L//o))^/^'^. These limit- 
ing behaviors have been inferred from asymptotic expansions of Eq. ()14p . 





FIG. 2: Dimensionless MF scaling function ^^^{x^^.uq = 0) = /c-L^(vo/90)^/2 (see Eq. (fTO]l ) 
with x^^^ = X = ro-L^ plotted together with the renormalized mean field scaling function 

f(.L^ = ^RMFf^^RMF^ (^ggg ^^^^ ^^^^-^ ^-^j^ ^RMF ^^^^^2 ^^^^^ ^v/Zq = 400. 

^MF^^MF ^ oo^o) ~ 8(x*^^)3/2e-2{^'^'^)'/' (thin dash-dotted line) and 7?*^^(x^^^ ^ 0,0) ~ 
2.76— 0.5x^^^ (thin dashed line). The asymptotic behavior oi^^^^^{x^^^^ , 0) can be obtained from 
the one of (x*^^ , 0) by multiplying the ordinate by the factor (28/(87rV3))^/^(ln(L//o))^/^; the 
abscissa remains the same. 
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FIG. 3: Experimental data from Ref. 4] for the scaling functions -d = fcL^ for the Casimir force in 
^He-^He films of thicknesses L along various paths of fixed ^He concentration (given in the figure) 
close to the tricritical concentration Xf = 0.672. The scaling variable is in units of A. The solid 
line corresponds to the tricritical mean-field scaling function [4] calculated for rg = (i.e., a = 
in Eq. (jlSp ) and suitably adjusted (see the main text); t = (T — Tt)/Tt. 
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FIG. 4: Bulk phase diagram for the VBEG model obtained whithin MFT for K/J = 0.5 and 
A^'-* / J = —3 exhibiting the line Ts{X) of continuous superfluid transitions in the bulk (long-dashed 
line), the phase separation curves (solid lines), the tricritical point A = {Tt/Ts{0) = 2/3, Xt = 1/3). 
In a semi-infinite system there is a (short-dashed) line of continuous surface transitions which 
merges with the line Ts{X) of bulk critical points at the special transition point S = (Ts/Ts{0) ~ 
0.759,^5 ~ 0.241). Upon crossing this surface transition line a thin film near the surface becomes 
superfluid although the bulk remains a normal fluid. Vertical lines represent thermodynamic paths 
along which the Casimir force has been calculated (see, c.f.. Fig. [5]). • (A) state points which 
will be considered in Fig. [71 
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y = tL 

FIG. 5: Dimensionless scaling function = tL) = fcL^, with t = {T — Tt) /Tt and Z = 20 for the 
Casimir force calculated within MFT for the VBEG model along the paths of fixed concentration 
of '^He shown in Fig. [H Dots indicate the corresponding onset temperature Ts{X) of superfluidity 
at the line of bulk critical points. The full line for y < corresponds to the temperatures of the 
onset of the first-order phase separation in the bulk (see Fig. U]). In view of, c.f., Fig. [9] we note 
that the curves might still shift if calculated for larger values of L. 
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FIG. 6: Bulk phase diagram for the VBEG model in the (A, T) plane obtained whithin MFT 
for the same set of parameters as in Fig. 21 The long-dashed coexistence line corresponds to the 
continuous superfluid transitions whereas the solid coexistence line corresponds to the curve of 
first-order phase separation. As indicated in the inset g and t are the two relevant scaling variables 
(compare Eq. (jl7p ): the line g' = is tangential to the coexistence line at the tricritical point where 
the lines of first- and second-order transitions merge. Note that according to Eq. (jl7p along the 
line 5 = one has (A — At)/{kBTt) = —a't and along the line t = one has g = {A — At) / {kBTt). 
Three thermodynamic paths of constant concentration are shown: X = Xf, X = X^ — 0.005 (upper 
line), and X = Xt + 0.005 (lower line). We note that along the paths of constant concentration 
both scaling variable t and g vary; however, the variation of t is more pronounced so that within 
a rough approximation g can be considered to be constant along each path. 
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FIG. 7: (a) '^He concentration profile X{1) = 1 — Qi and (b) superfluid OP profile mi for a VBEG 
film of thickness Z = 60 for = 0.5 J, A*^')/J = —3, and A^^'^ /J = At/ J ~ 0.61 corresponding to 
the state points •, and ▲ indicated in Fig. UJ t = (T — Tt)/Tt. 
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FIG. 8: (a) '^He concentration profile X(l) = 1 — Qi and (b) superfluid OP profile mi for a VBEG 
film of width Z = 60 for K = 0.5 J, A^'-'/J = — oo, and A^^^/J = +oo corresponding to the state 
points •, and ▲ indicated in Fig. d] with t = {T — Tt)/Tt. 
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FIG. 9: Mean-field scaling function -i^oCy = r(L/^(]")^) = fcL'^ for the limiting case of the VBEG 
model (symbols) corresponding to pure ^He and various film thicknesses L with t = (T — T\)/T\. 
The full curve corresponds to the scaling function 'dQ^{y) obtained from the continuum 0(2) LG 
theory within MFT (Eqs. (j55|) and (j56|) ) whith the amplitude = Am{L) and the position of 
the minimum Um = ym{L) determined in such a way as to provide the best fit to "i^o from the 
VBEG model; for further details see the main text. With this rescaling the continuum theory 
provides a very good fit (here shown only for L = 60) to the numerical data. The insets show 
the L-dependence of and Um used as fitting parameters. The dashed line in the inset for 
ym{L) indicates the limiting value ym = — vr^ predicted by th LG model. Surprisingly, scaling - 
corresponding to Z-independent Am and ym - is not yet attained by the numerical data of the 
VBEG model even for thick slabs with L ~ 60. 



41 



0.2 0.4 0.6 0.8 

z/L 



FIG. 10: Mean-field OP profiles (normalized to the corresponding bulk values nifj) across slabs of 
thickness L calculated from the limiting case of the VBEG model (symbols, L = 150) and from the 



continuum 0(2) LG theory (lines, see Eqs. (202) and (203) in Ref. [3^) for a selection of the scaling 
variable y = t{L/£^qY/'^ below the shifted critical point of the film (corresponding to y = ym = — tt^, 
see the main text). For y sufficiently negative m{z ^ a) — rrih ~ exp(— 2:/^(r < 0)) in the middle 
of the slab. This allows one to infer = ^(r < 0)(-r)i/2 ~ 0.29 so that ^+ = ^2^^ ~ 0.41. 
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FIG. 11: The adjusted scaling function d{y) (see the main text) for the VBEG model within 
MFT compared with the corresponding experimental curve obtained along the path of fixed 
tricritical concentration X = Xt ^ 0.672 of ^He. ^{y) is obtained from in Fig. 5 by rescaling 
the amplitudes of -d and y such that there is agreement between the experimental data for X = 0.672 
at y = and with respect to the positions of the maximum. The VBEG curve for X = Xf — 0.01 
agrees with the experimental data for nominally X = Xt even better. Both theoretical curves 
coincide for y < 0. 
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FIG. 12: Normalized mean-field scaling function i9o(y) for the limiting case of the VBEG model 



(on a lattice with L = 60) corresponding to pure '^He compared with the experimental data {'&) 
2] in terms of the proper scaling variable y = t(^L / ^q)^/'^ using (^Q')exp = 

1.43A for pure He 
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and u = 0.67. These are the universal forms of the scaling function i?o- The inset on the left 
shows a magnification of the main plot close to the minimum. According to the analysis presented 
in Subsec. IIVBI (see also Fig. [9|) the position ym{L) of the minimum of the theoretical curve in 
the scaling limit L = oo approaches the value — vr^. In the inset on the right the experimental 
data (diamonds) above the critical temperature are compared with the scaling functions for the 
three-dimensional XY model in a slab obtained from the e-expansion (see the main text: the 
solid (dashed) line corresponds to the [1,0] ([0,1]) Fade approximant). Due to the experimental 
resolution {'&)exp takes only discretized values. 
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